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Abstract. A hybrid lattice-statistical model of doubly decorated two-dimensional lattices, 
which have localized Ising spins at its nodal sites and itinerant electrons delocalized 
over decorating sites, is exactly solved with the help of a generalized decoration-iteration 
transformation. Under the assumption of a quarter filling of each couple of the decorating 
sites, the ground state constitutes either spontaneously long-range ordered ferromagnetic 
or ferrimagnetic phase in dependence on whether the ferromagnetic or antiferromagnetic 
interaction between the localized Ising spins and itinerant electrons is considered. The critical 
temperature of the spontaneously long-range ordered phases monotonically increases upon 
strengthening the ratio between the kinetic term and the Ising-type exchange interaction. 



1. Introduction 

Exactly soluble lattice-statistical models traditionally attract appreciable scientific interest as 
they offer a valuable insight into diverse aspects of cooperative phenomena p]. It should 
be mentioned, however, that sophisticated mathematical methods must be usually employed 
when searching for an exact treatment of even relatively simple interacting many-body systems 
and consequently, a theoretical treatment of more realistic or more complex models is often 
accompanied with a substantial increase of computational difficulties. The mapping technique 
based on generalised algebraic transformations belongs to the simplest mathematical methods, 
which allow to obtain the exact solution of a more complicated model from a precise mapping 
relationship with a simpler exactly solved model [2, 3j. Recently, this approach has been applied 
to an intriguing diamond chain model of interacting spin-electron system [HE]. The purpose of 
this work is to treat exactly similar two-dimensional (2D) model, which should provide a deeper 
insight into phase transitions and critical phenomena of interacting spin-electron systems. 



2. Model and its exact solution 

Let us consider a hybrid lattice-statistical model of interacting spin-electron system on doubly 
decorated 2D lattices, which have one localized Ising spin at each nodal site and one itinerant 
electron delocalized over each couple of decorating sites. The magnetic structure of the model 
under investigation is schematically illustrated in figure [T] on the particular example of the 
doubly decorated square lattice. For further convenience, let us define the total Hamiltonian as 
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Figure 1. A part of the doubly decorated 
square lattice, which has at each nodal site 
(full circle) one localized Ising spin and 
at each couple of decorating sites (empty 
circles) one itinerant electron. The ellipse 
demarcates unit cell described by the kth. 
bond Hamiltonian given by Eq. ([T]). 



a sum over bond Hamiltonians Ti. = X^fc^fc) where each bond Hamiltonian Ti./. involves all the 
interaction terms of one itinerant electron from the kth. couple of the decorating sites 



t 



J 

2 



(4. 



(1) 



In above, ^ and c^^^ (a = 1,2, 7 =T,i) denote usual creation and annihilation fermionic 



"ka,"/ 



operators and (t|^ is the standard spin-1/2 operator with the eigenvalues cr|^ = ±1/2. The 
transfer integral t takes into account a kinetic energy of a single electron delocalized over a 
couple of the decorating sites and the exchange integral J describes the Ising-type interaction 
between the delocalized electron and its two nearest-neighbouring localized Ising spins. 

The crucial step of our calculation represents an evaluation of the partition function. A 
validity of the commutation relation between different bond Hamiltonians ['Hi,Ti.j] = allows 
a straightforward factorization of the total partition function Z into a product performed over 
bond partition functions Z^. 
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where /? = l/^k^T), k^ is Boltzmann's constant, T is the absolute temperature, is the total 
number of the Ising spins (i.e. the nodal lattice sites) and q is their coordination number (i.e. the 
number of nearest neighbours). Next, the symbols Trj^.j and Trj^.} denote a trace over degrees 
of freedom of all Ising spins and itinerant electrons, respectively, while the symbol T^^Cki,Ck2 stands 
for a trace over degrees of freedom of the itinerant electron from the kth. couple of decorating 
sites. After the elementary diagonalisation of the bond Hamiltonian 7ik, one arrives at the 
explicit expression of the bond partition function Z^, which can be eventually replaced through 
the appropriately chosen generalised decoration-iteration transformation [21 [3] 
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The 'self-consistency' condition of the algebraic transformation ^ requires that this mapping 
relationship must hold independently of the spin states of two Ising spins a^^ and involved 
therein, which directly determines so far not specified transformation parameters A and R as 
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At this stage, a substitution of the generahsed decoration-iteration transformation ([3]) into 
Eq. ([2]) yields in turn a simple mapping relation between the partition function of the interacting 
spin-electron system on the doubly decorated 2D lattice and the partition function of the simple 
spin- 1/2 Ising model on the corresponding undecorated lattice 

Z{P,J,t)=A^''/'ZiM{f3,R). (5) 

It can be easily understood from Eqs. dH) and ^ that the mapping parameter A cannot cause 
a non-analytic behaviour of the partition function Z and hence, the investigated spin-electron 
system becomes critical if and only if the corresponding spin-1/2 Ising model becomes critical as 
well. Accordingly, the lines of critical points can readily be obtained from a comparison of the 
effective temperature-dependent coupling PR with the relevant critical points of the spin-1/2 
Ising model on the corresponding undecorated lattices that are given by pcR = 21n(2 + ^/3), 
21n(l -|- \/2), and In 3 for the honeycomb, square, and triangular lattices [6], respectively. 

Other thermodynamic quantities can be now easily derived from the mapping relation ([5]) 
between the partition functions. For instance, both sublattice magnetisations can be calculated 
by combining the mapping relation ([5|) with the exact Callen-Suzuki identity. As a result, the 
sublattice magnetisations rrii and rrif. per one Ising spin and per one itinerant electron read 

nii = miM(/3-R), rUe = tanh (/3J/4) niiMiPR), (6) 

where miMi(3R) denotes the spontaneous magnetisation of the spin-1/2 Ising model on the 
corresponding undecorated lattice that is known for several planar lattices [6]. 

3. Results and discussion 

Let us proceed to a discussion of the most interesting results obtained in the preceding section. 
First, it should be realized that the effective coupling ^ of the spin-1/2 Ising model on the 
corresponding undecorated lattice is always positive, i.e. (3R > 0, which means that the 
interacting spin-electron system is effectively mapped to the ferromagnetic spin-1/2 Ising model. 
Hence, it directly follows from Eq. ([6]) that the ground state constitutes either the classical 
ferromagnetic phase with fully saturated and identically oriented sublattice magnetisations 
rui = nie = 1/2 on assumption that J > 0, or the classical ferrimagnetic phase with fully 
saturated sublattice magnetisations oriented opposite one to each other rui = —me = 1/2 
if J < 0. Both the spontaneously long-range ordered phases exhibit completely the same 
critical behaviour as a result of the invariance of the effective coupling PR with respect to 
the transformation J ^ — J. 

The reduced critical temperature of the spontaneously long-range ordered phases is plotted 
in Fig. [2] against the ratio t/\J\ between the hopping term and the exchange integral. As one 
can see, the critical temperature rises steadily when increasing a relative strength of the kinetic 
term until it asymptotically tends towards its maximum value achieved in the limit t/\J\ — > oo. 
In the limit t/\J\ 0, the observed zero critical temperature is consistent with a non-magnetic 
character of one of two decorating sites on each bond of the doubly decorated lattice. On 
the other hand, the highest asymptotic values of the critical temperatures for the interacting 
spin-electron system on the doubly decorated honeycomb, square and triangular lattices 

kBTc 1 

'^^'^ 41n(2 + V3 + \/6 + 4^/3)' ^^^^ 41n(l + ^/2 + \/2 + 2^/2) ' '^^'^ 41n(y2 + ^/3)' 

are exactly a half of the relevant critical temperatures of the spin-1/2 Ising model on singly 
decorated honeycomb, square, and triangular lattices, respectively. This observation would 
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Figure 2. The dimensionless critical tempera- 
ture as a function of the relative strength of the 
kinetic term for doubly decorated honeycomb, 
square, and triangular lattices. 
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Figure 3. Temperature variations of the 
total and both sublattice magnetisations for 
the doubly decorated square lattice at three 
different values of the ratio t/ J . 



suggest that a delocalization of itinerant electrons generally reduces the critical temperature, 
because the localized Ising spins effectively feel due to the hopping process less than a half of 
the total magnetic moment of each itinerant electron. 

Thermal variations of the total and sublattice magnetisations are displayed in Fig. [3] for three 
different values of a relative strength of the kinetic term. It is quite obvious from this figure 
that both sublattice magnetisations exhibit quite similar temperature dependencies when the 
sublattice magnetization of the localized Ising spins lies just slightly above the sublattice 
magnetisation rUe of the itinerant electrons. A closer analysis reveals that both sublattice 
magnetisations tend to zero in a vicinity of the critical temperature with the critical exponent 
from the standard Ising universality class. 

To conclude, we have found the exact solution for the hybrid model of the interacting spin- 
electron system on doubly decorated 2D lattices. Under the assumption of quarter filling of 
each couple of the decorating sites, the ground state constitutes either spontaneously long- 
range ordered ferromagnetic or ferrimagnetic phase in dependence on whether the ferromagnetic 
or antiferromagnetic interaction between the localized Ising spins and itinerant electrons is 
assumed. It has been shown that the critical temperature of spontaneously long-range ordered 
phases monotonically increases upon strengthening the ratio between the kinetic term and the 
exchange interaction. The work on an analogous 2D hybrid model of interacting spin-electron 
system with two electrons per each couple of the decorating sites is in progress [7j. 
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